Cooperative effect of load and disorder in thermally 
activated rupture of a 2d random fuse network 



A. Guarino^, L. VaneP, R. Scorretti^ and S. Ciliberto^ 

^ Universite de la Polyncsic Frangaise, BP 6570 Faaa Aeroport, Tahiti, French 
Polynesia 

^ Laboratoire de physique, CNRS UMR 5672, Ecole Normale Superieure de Lyon, 46 
allee d'ltalie, 69364 Lyon Cedex 07, France 

3 CEGELY - UMR 5005, Universite Lyon, 1 43, bid du 11 Nov. 1918, 69622 
Villeurbanne Cedex, France 

E-mail: Loic.VanelOens-lyon.fr 

Abstract. A random fuse network, or equivalently a 2d spring network with 
quenched disorder, is submitted to a constant load and thermal noise, and studied 
by numerical simulations. Rupture is thermally activated and the lifetime follows 
an Arrhenius law where the energy barrier is reduced by disorder. Due to the 
non-homogenous distribution of forces from stress concentration at microcracks' tips, 
spatial correlations between rupture events appear, but they do not affect the energy 
barrier's dependence on disorder, they affect only the coupling between disorder and 
the applied load. 

PACS numbers: 61.43.-j, 05.10.-a, 05.70.Ln, 62.20.Mk 

The effect of quenclied disorder on dynamics is a recurring problem in many pfiysical 
systems witli elastic interactions. The motion of vortex lines in supraconductors, charge- 
density waves in Bragg glasses, magnetic domains walls, or contact lines in wetting 
show a competition between elastic interactions and pinning by disorder |2] . While 
many studies have focused on systems driven above a critical depininng threshold, an 
important issue remains to understand the sub-critical regime, when thermally activated 
creep motion occurs [3 E] . 

Rupture in disordered brittle solids falls in the same class of problems jS]. Elastic 
interactions tend to make a crack propagate in a straight direction while disorder creates 
roughness or causes spatially diffuse damage [71 |H]. In the sub-critical rupture regime, 
a very important quantity for safety reasons is the lifetime, i.e. the mean time for a 
sample to break under a prescribed load. The lifetime follows an Arrehnius law llUj. 
but thermal noise is generally too small compared to recent theoretical estimates of 
the energy barrier ^T] ^] to explain experimental observations in heterogenous 
materials [HI [13 HEj- For athermal systems, disorder actually reduces the energy 
barrier and can be seen as an effective temperature j^. In order to clarify the role 
of disorder in thermal systems, one- dimensional Thermal and Disordered Fiber Bundles 
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Models (l(i-TDFBM) have been introduced to model the thermally activated rupture 
of an heterogenous material submitted to a constant external load [THl IT^ 1^ 1^ 1^. 
The TDFBM considers an elastic system in equilibrium at constant temperature where 
statistical force fluctuations occur in time due to thermal noise. This is very different 
from previous thermal random fuse models j2S! where rupture results from an increase 
in fuse temperature due to dissipation through a generalized Joule effect until reaching a 
critical melting temperature. In a TDFBM, elastic energy is the equivalent of dissipation 
in the thermal random fuse model but does not cause rupture when the system is at 
mechanical equilibrium; rupture is caused by elastic force fluctuations analogous to 
Nyquist noise. 

One problem with the Irf-TDFBM investigated up to now is that the load is shared 
equally among all the unbroken fibers. This is not a realistic load-sharing rule for 
experimental geometries where elasticity cause stress concentration at microcracks' tips 
and lead to a non-uniform redistribution of stress. In this letter, we show that spatial 
correlations between rupture events in 2d do not affect the dependence of the energy 
barrier on disorder but only the coupling between disorder and the applied load. 

First, we discuss briefly results obtained by several authors [IHl CHI EDI EH 122] on the 
l(i-TDFBM. The system considered is made of a set of parallel fibers, each carrying 
an initial force /o and behaving as a linear elastic spring with unity stiffness. Each fiber 
j can carry a maximum force fc^ before it breaks. Quenched disorder is introduced 
in the system by distributing thresholds fc^ according to a gaussian distribution of 
mean < /i"'^ >= 1 and variance T^; for each fiber, the value /i"''' is a time-independent 
constant. Contrary to the case of non-thermal l(i-DFBM where the system evolves due 
to a progressive increase in total current, we consider that the total force applied to 
the l(i-TDFBM is kept constant. Dynamics is introduced in the system by introducing 
fluctuations in spring forces due to thermal noise. We write fj the average force on fiber 
j. The fluctuations in force 6fj that occur in time on fiber j are assumed to follow a 
gaussian probability distribution with mean value and variance T, where T represents 
the thermodynamical temperature in unit of square force. When the total force on a 
fiber fj + 6fj is larger than the threshold fc''\ the fiber breaks. The remaining fibers 
share equally the total force: this is a so-called democratic model. The bundle will break 
completely as soon as the average force on each fiber exceeds the breaking threshold. 
Roux has shown that the mean time to break the first fiber follows an Arrhenius law 
where disorder acts as an additive temperature [T^ : 



with /o the initial force carried by each fiber of the bundle. In the general case where 
many fibers break before total rupture of the bundle, the lifetime obeys an Arrhenius 
law with a different general form : 




(1) 




(2) 



Cooperative effect of load and disorder in thermally activated rupture 



3 



An approximate expression of U for low disorder is jl9[ 1201 ■ 

t/ = ^(l-/o-/5iV^)' (3) 

with /9i = A/vr/2. For higher disorder, Politi et al pT^ have shown numerically that U is 
determined with a very good approximation by the minimum value of the rate at which 
fibers break. More precisely, if n is the number of broken fibers and N the total number 
of fibers then ^ = n/N is the fraction of broken fibers and $ its time derivative. The 
value $* for which $ is minimum obeys an implicit equation pT| |22] : 

exp [irfc^(2<l>*)] (1 - = a 

where irfc is the inverse function of the complementary error function and a = 
fo/y/2nTd. Then, U can be approximated as j21]: 

f/ = ^(l-/o-/?2V^)' With /?2(«) = «v^Y^ + v^irfc(2$*) (4) 

Note that the coupling coefficient P2 depends on both disorder and applied load 
/o. Eq.(jD) predicts that the variation of U with Td for a fixed /o is non monotonous 
[22] • However, when a > 1/4 (this condition corresponds to $* < 1/2 which would 
be expected in practice for most materials), U decreases when Td is increased with 
/o constant. In that case, the function /32 has a lower value /32(l/4) = /3i/2 and an 
asymptotic behavior for large a ($* 0), (32{<^) — v21na + 1 / V2 In a [2X. 

To study the effect of a non-uniform force redistribution on the rupture dynamics of 
the TDFBM, one could keep a Id geometry and introduce a finite range of interaction 
between fibers |25j. Instead, we work directly in a 2d geometry more closely related 
to a real experiment. The above described Id-TDFBM is equivalent to a system of 
parallel fuses where forces are transformed in currents and displacements in electric 
potentials. A 2d square fuse network is then equivalent to a square lattice of springs in 
antiplane deformation. Specifically, each node of the N x N nodes square lattice can 
move along an axis perpendicular to the plane of springs at rest. A constant force F 
is applied at two opposite sides of the lattice in antiplane configuration. In the initial 
equilibrium configuration of the lattice, the springs submitted to a load /o = F/N are 
called "parallel" springs while the unloaded springs (zero force) are called "orthogonal" 
springs. 

Like in Id, the fluctuations in force 6fj on spring j follow a gaussian probability 
distribution with mean value and variance T, and the rupture thresholds /i"''' follow 
a gaussian distribution of mean < fc ^ >= 1 and variance Td- The time scale in the 
simulation is the (constant) time between two configurations of force fluctuations in 
the system. The whole network is a square with sides 100 springs wide; thus, the 
lattice contains about 2. 10'^ springs. Whenever we choose rupture thresholds from the 
gaussian distribution, there is a non-zero probability to obtain a negative threshold : 
P<o = |erfc (< fc > / V2Td) . For a system with 210^ springs, we can safely consider 
that no spring will spontaneously break when there is no load (/o = 0) and at zero 
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temperature if P<o < 10 ^, thus < 0.055 when < fc >= 1. In practice, we will work 
with Td < 0.05. 

First, let us consider a lattice with no disorder in the rupture thresholds and all 
the springs initially intact. When thermal noise is on, some springs will start to break. 
Due to stress concentration effects, as soon as one of the parallel springs is broken, 
the neighboring springs are submitted to a larger force. If that new force exceeds the 
rupture threshold then the neighboring spring will also break, and the force on the next 
neighbor will be even higher. This process will result in the rupture of the whole network 
in an avalanche started from a single rupture event. Numerically, this will occur in our 
lattice as soon as /o > , where (— 0.785) is the critical threshold of the homogenous 
lattice at T = when one parallel spring is broken. In that case, the rupture time will 
be directly related to the probability of breaking a single spring in the network and, if 
the lattice is disordered, we will recover essentially the result given by eq. ((H) [THj. In 
the rest of the paper, we will be interested only in the case where several springs break 
before the final avalanche occurs, i.e. when /o < f^. 
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Figure 1. Scaled rupture time versus Arrlienius factor U/T. Various symbols 
correspond to different values of disorder and applied force. The solid line corresponds 
to logy = X. 

For a fixed value of the applied force /o and a fixed value of the disorder T^, we 
measure as a function of temperature T the lifetime, i.e. the rupture time averaged over 
as many as 100 numerical experiments. We find that lifetime follows an Arrhenius law 
and obtain numerically f/(/o, T^) as defined by eq ©. We see on figure[T]the collapse of 
the data for a range of values 0.04 < fo < 0.77 and 10~^ < < 0.05. Some of the data 
points are more scattered around the expected behavior (solid line) than others because 
the ratio of the standard deviation over the mean of the rupture time increases when T 
becomes small. This property, already mentioned in jTH], makes numerical convergence 
of the mean difficult in some cases. 

To compare the barrier U{fo, Tf) for the 2d geometry with the one of the Id- 
TDFBM, we plot on figure v^/(l - fo) as a function of /?2V^/(1 - fo)- We 
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Figure 2. Comparison between the numerical value of U in 2d and eq. in Id (solid 
line). Inset : same scaling, but using the numerical value Uq = [/(/o,0) for the 2d 
lattice. 




Figure 3. Collapse of all the numerical data using an effective value Pe which is only 
a frmction of /o but not of T^, contrary to the case of f32 which is a function of T^. 

clearly see that the functional form for /?2 is not correct, but also that the data for 
various values of disorder do not rescale very well. The first immediate reason for the 
discrepancy is that for = 0, the energy barrier of the 2d network Uq = U{fo,0) 
(determined numerically) is not (1 — /o)^/2 as in Id. This is due to the preferential 
redistribution on the nearest neighbors of the force carried by a fiber before rupture. 
Taking into account the effective energy barrier in 2d does not improve the comparison 
with the Id-TDFBM. After replacing 1 — /q by \/2Uq in eq. (jH), we see in the inset of 
figure 121 that it is not enough to get a good collapse of all the data on the theoretical 
prediction (solid line). 

We find that the barrier decreases with disorder following the linear curve : 
A/f/(/o, Tf) = y/Uo — by/Id for 10~^ < Td < 0.05 and a fixed value of /q. Thus, it 
turns out that eq. (jH)) is a much better functional form than eq. (j3)), even though it is 
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Figure 4. Numerical values of f3e as a function of fo (circles) . The dashed and dotted 
lines correspond to f32 as predicted by eq. when Td — 0.001 and Td = 0.01. Close 
to f^, the values found for Pe are close to values predicted by eq. Q|. 



an approximate expression in Id. Looking at eq. Q or eq. we can make an analogy 
with the Id case and say that the second coefficient b corresponds to an effective value 
Pe/V2 which is now only a function of /q. On figure El we see the collapse of all the 
data when we plot ^^U/Uq as a function of l3e{fo)\/Td/ V^Uq. 

The coupling coefficient Pe increases almost linearly with /o up to values close to 
f^ (figure Ej). However, when fo gets very close to f^, there is an abrupt decrease in the 
value of /3e- This is related to the fact that rupture is now controlled by a single event 
as in eq. (|T)). Indeed, although eq. (H)) does not follow the general scaling property of 
eq.((21), we can estimate a value Pe for each temperature value used in the simulation. 
The average value found for (5^ from eq. (Q) (square symbols in figured with an error bar 
corresponding to variations with temperature) is a reasonable estimate of the numerical 
value. 

The functional behavior of (3e is very different from that of the Id model (eq. (^) 
where P2 depends on both fo and T^. As an example, we plot on figure |ll/32(/o) for fixed 
values of T^. Not only the functional dependence is clearly different from the numerical 
estimate Pe{fo) but also P2 decreases with at fixed fo. In that sense, the load and 
the disorder do not act cooperatively in Id. 

The key point in 2d is that the spatial correlations between rupture events depend 
on the strength of stress intensification. This is illustrated on figure El a) to c) which 
shows the broken fibers just before the final avalanche for different values of fo and a 
fixed value of Td. For very small loads, damage is scattered everywhere in the sample. 
At higher loads, damage becomes less scattered and growth of straight cracks occurs. 
Finally, figure Eli) shows that for a load close to the critical threshold, only very few 
events occur. A similar transition was observed for zero disorder or annealed disorder 
models with power law rate of rupture |2nil2Zll2Hl- In contrast to these models where 
the transition occurs by changing the exponent of the power law, we observe a transition 
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Figure 5. Image of the broken fibers (black dots) just before the final rupture for 
Td = 0.01 and /o = 0.05 (top left), /o = 0.1 (top right), fo = 0.3 (bottom left); 
Td — 0.004 and fo = 0.75 (bottom right). The arrow shows the loading direction. 

resulting from the competition between stress intensification and quenched disorder. 

To understand this transition in our model, let us consider the increase in force 
due to stress intensification when a spring breaks. If the increase is small compared 
to y/Td, there will be very little spatial correlation between rupture events occurring 
preferentially at the weakest springs. For a given disorder, there is always a force fo small 
enough to observe this rupture regime similar to the Irf-TDFBM case. On the contrary, 
if the increase due to stress intensification is large compared to y/Td, it is easier to 
break a spring next to an already broken spring, and the rupture will proceed mainly by 
growth of multiple cracks. In spite of very different regimes of spatial correlation between 
rupture events, we have the remarkable result that the energy barrier dependence on 
Td is unchanged. Spatial correlations only affects the coupling coefficient jSg, increasing 
quasi- linearly with fo and independent of T^. 

The multiplicative amplification of disorder due to (3e is a mechanism that will 
create a load-dependent reduction of the energy barrier in thermally activated rupture. 
It will have an effect on the order of magnitude and load-dependence of the rupture 
time which could help understanding experiments in heterogeneous materials [Tl] . 

In conclusion, we have studied thermally activated rupture of a 2d elastic spring 
network submitted to a constant load and thermal noise. We find that spatial 
correlations between rupture events are controlled by a competition between quenched 
disorder and force inhomogeneities due to stress concentration. For low spatial 
correlations, the energy barrier scales naturally like in the Id model. Remarkably, 
the appearance of spatial correlations does not affect the functional dependence of the 
energy barrier on disorder, but only the coupling coefficient f3e which is independent 
of disorder and increases quasi-linearly with the applied load fo- This is an important 
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result showing that the apphed load contribute to amplify in a cooperative way the 
effect of disorder on the lifetime. The observed cooperative effects of load and disorder 
in 2d subcritical rupture could be relevant to the creep regime of other physical systems 
with elastic interactions ^1 12] and also to crackling noise [22] • 
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